Search

We present the integrated 3-point correlation functions (3PCF) involving both the cosmic shear and the galaxy density fields. These are a set of higher-order statistics that describe the modulation of local 2-point correlation functions (2PCF) by large-scale features in the fields, and which are easy to measure from galaxy imaging surveys. Based on previous works on the shear-only integrated 3PCF, we develop the theoretical framework for modelling 5 new statistics involving the galaxy field and its cross-correlations with cosmic shear. Using realistic galaxy and cosmic shear mocks from simulations, we determine the regime of validity of our models based on leading-order standard perturbation theory with an MCMC analysis that recovers unbiased constraints of the amplitude of fluctuations parameter $A_s$ and the linear and quadratic galaxy bias parameters $b_1$ and $b_2$. Using Fisher matrix forecasts for a DES-Y3-like survey, relative to baseline analyses with conventional 3$\times$2PCFs, we find that the addition of the shear-only integrated 3PCF can improve cosmological parameter constraints by $20-40\%$. The subsequent addition of the new statistics introduced in this paper can lead to further improvements of $10-20\%$, even when utilizing only conservatively large scales where the tree-level models are valid. Our results motivate future work on the galaxy and shear integrated 3PCFs, which offer a practical way to extend standard analyses based on 3$\times$2PCFs to systematically probe the non-Gaussian information content of cosmic density fields., Comment: 19 pages, 8 figures + appendix. Comments are welcome!

We build a deep learning framework that connects the local formation process of dark matter halos to the halo bias. We train a convolutional neural network (CNN) to predict the final mass and concentration of dark matter halos from the initial conditions. The CNN is then used as a surrogate model to derive the response of the halos' mass and concentration to long-wavelength perturbations in the initial conditions, and consequently the halo bias parameters following the "response bias" definition. The CNN correctly predicts how the local properties of dark matter halos respond to changes in the large-scale environment, despite no explicit knowledge of halo bias being provided during training. We show that the CNN recovers the known trends for the linear and second-order density bias parameters $b_1$ and $b_2$, as well as for the local primordial non-Gaussianity linear bias parameter $b_\phi$. The expected secondary assembly bias dependence on halo concentration is also recovered by the CNN: at fixed mass, halo concentration has only a mild impact on $b_1$, but a strong impact on $b_\phi$. Our framework opens a new window for discovering which physical aspects of the halo's Lagrangian patch determine assembly bias, which in turn can inform physical models of halo formation and bias., Comment: 11 pages, 5 figures, to be submitted to MNRAS, comments welcome

The integrated shear 3-point correlation function $\zeta_{\pm}$ measures the correlation between the local shear 2-point function $\xi_{\pm}$ and the 1-point shear aperture mass in patches of the sky. Unlike other higher-order statistics, $\zeta_{\pm}$ can be efficiently measured from cosmic shear data, and it admits accurate theory predictions on a wide range of scales as a function of cosmological and baryonic feedback parameters. Here, we develop and test a likelihood analysis pipeline for cosmological constraints using $\zeta_{\pm}$. We incorporate treatment of systematic effects from photometric redshift uncertainties, shear calibration bias and galaxy intrinsic alignments. We also develop an accurate neural-network emulator for fast theory predictions in MCMC parameter inference analyses. We test our pipeline using realistic cosmic shear maps based on $N$-body simulations with a DES Y3-like footprint, mask and source tomographic bins, finding unbiased parameter constraints. Relative to $\xi_{\pm}$-only, adding $\zeta_{\pm}$ can lead to $\approx 10-25\%$ improvements on the constraints of parameters like $A_s$ (or $\sigma_8$) and $w_0$. We find no evidence in $\xi_{\pm} + \zeta_{\pm}$ constraints of a significant mitigation of the impact of systematics. We also investigate the impact of the size of the apertures where $\zeta_{\pm}$ is measured, and of the strategy to estimate the covariance matrix ($N$-body vs. lognormal). Our analysis solidifies the strong potential of the $\zeta_{\pm}$ statistic and puts forward a pipeline that can be readily used to improve cosmological constraints using real cosmic shear data., Comment: 21 pages, 11 figures, 3 tables. Comments welcome

We discuss the potential of the multi-tracer technique to improve observational constraints of the local primordial non-Gaussianity (PNG) parameter $f_{\rm NL}$ from the galaxy power spectrum. For two galaxy samples $A$ and $B$, the constraining power is $\propto |b_1^B b_\phi^A - b_1^Ab_\phi^B|$, where $b_1$ and $b_\phi$ are the linear and PNG galaxy bias parameters. We show this allows for significantly improved constraints compared to the traditional expectation $\propto |b_1^A - b_1^B|$ based on naive universality-like relations where $b_\phi \propto b_1$. Using IllustrisTNG galaxy simulation data, we find that different equal galaxy number splits of the full sample lead to different $|b_1^B b_\phi^A - b_1^Ab_\phi^B|$, and thus have different constraining power. Of all of the strategies explored, splitting by $g-r$ color is the most promising, more than doubling the significance of detecting $f_{\rm NL}b_\phi \neq 0$. Importantly, since these are constraints on $f_{\rm NL}b_\phi$ and not $f_{\rm NL}$, they do not require priors on the $b_\phi(b_1)$ relation. For direct constraints on $f_{\rm NL}$, we show that multi-tracer constraints can be significantly more robust than single-tracer to $b_\phi$ misspecifications and uncertainties; this relaxes the precision and accuracy requirements for $b_\phi$ priors. Our results present new opportunities to improve our chances to detect and robustly constrain $f_{\rm NL}$, and strongly motivate galaxy formation simulation campaigns to calibrate the $b_\phi(b_1)$ relation., Comment: 9 pages, 4 figures, 1 table. Comments welcomed! v2 matches version published in JCAP

We use the BOSS DR12 galaxy power spectrum to constrain compensated isocurvature perturbations (CIP), which are opposite-sign primordial baryon and dark matter perturbations that leave the total matter density unchanged. Long-wavelength CIP $\sigma(\vec{x})$ enter the galaxy density contrast as $\delta_g(\vec{x}) \supset b_\sigma\sigma(\vec{x})$, with $b_\sigma$ the linear CIP galaxy bias parameter. We parameterize the CIP spectra as $P_{\sigma\sigma} = A^2P_{\mathcal{R}\mathcal{R}}$ and $P_{\sigma\mathcal{R}} = \xi\sqrt{P_{\sigma\sigma}P_{\mathcal{R}\mathcal{R}}}$, where $A$ is the CIP amplitude and $\xi$ is the correlation with the curvature perturbations $\mathcal{R}$. We find a significance of detection of $Ab_\sigma \neq 0$ of $1.8\sigma$ for correlated ($\xi = 1$) and $3.7\sigma$ for uncorrelated ($\xi = 0$) CIP. Large-scale data systematics have a bigger impact for uncorrelated CIP, which may explain the large significance of detection. The constraints on $A$ depend on the assumed priors for the $b_\sigma$ parameter, which we estimate using separate universe simulations. Assuming $b_\sigma$ values representative of all halos we find $\sigma_A = 145$ for correlated CIP and $\sigma_{|A|} = 475$ for uncorrelated CIP. Our strongest uncorrelated CIP constraint is for $b_\sigma$ representative of the $33\%$ most concentrated halos, $\sigma_{|A|} = 197$, which is better than the current CMB bounds $|A| \lesssim 360$. We also discuss the impact of the local primordial non-Gaussianity parameter $f_{\rm NL}$ in CIP constraints. Our results demonstrate the power of galaxy data to place tight constraints on CIP, and motivate works to understand better the impact of data systematics, as well as to determine theory priors for $b_\sigma$., Comment: 11 pages, 4 figures, 3 tables. Comments welcomed!

We study the effect of super-sample covariance (SSC) on the power spectrum and higher-order statistics: bispectrum, halo mass function, and void size function. We also investigate the effect of SSC on the cross covariance between the statistics. We consider both the matter and halo fields. Higher-order statistics of the large-scale structure contain additional cosmological information beyond the power spectrum and are a powerful tool to constrain cosmology. They are a promising probe for ongoing and upcoming high precision cosmological surveys such as DESI, PFS, Rubin Observatory LSST, Euclid, SPHEREx, SKA, and Roman Space Telescope. Cosmological simulations used in modeling and validating these statistics often have sizes that are much smaller than the observed Universe. Density fluctuations on scales larger than the simulation box, known as super-sample modes, are not captured by the simulations and in turn can lead to inaccuracies in the covariance matrix. We compare the covariance measured using simulation boxes containing super-sample modes to those without. We also compare with the Separate Universe approach. We find that while the power spectrum, bispectrum and halo mass function show significant scale- or mass-dependent SSC, the void size function shows relatively small SSC. We also find significant SSC contributions to the cross covariances between the different statistics, implying that future joint-analyses will need to carefully take into consideration the effect of SSC. To enable further study of SSC, our simulations have been made publicly available at https://github.com/HalfDomeSims/ssc., Comment: 9+2 pages, 4 figures

We use $N$-body simulations to study halo assembly bias (i.e., the dependence of halo clustering on properties beyond total mass) in the density and primordial non-Gaussianity (PNG) linear bias parameters $b_1$ and $b_\phi$, respectively. We consider concentration, spin and sphericity as secondary halo properties, for which we find a clear detection of assembly bias for $b_1$ and $b_\phi$. At fixed total mass, halo spin and sphericity impact $b_1$ and $b_\phi$ in a similar manner, roughly preserving the shape of the linear $b_\phi(b_1)$ relation satisfied by the global halo population. Halo concentration, however, drives $b_1$ and $b_\phi$ in opposite directions. This induces significant changes to the $b_\phi(b_1)$ relation, with higher concentration halos having higher amplitude of $b_\phi(b_1)$. For $z=0.5$ and $b_1 \approx 2$ in particular, the population comprising either all halos, those with the $33\%$ lowest or those with the $33\%$ highest concentrations have a PNG bias of $b_\phi \approx 3$, $b_\phi \approx -1$ and $b_\phi \approx 9$, respectively. Varying the halo concentration can make $b_\phi$ very small and even change its sign. These results have important ramifications for galaxy clustering constraints of the local PNG parameter $f_{\rm NL}$ that assume fixed forms for the $b_\phi(b_1)$ relation. We illustrate the significant impact of halo assembly bias in actual data using the BOSS DR12 galaxy power spectrum: assuming that BOSS galaxies are representative of all halos, the $33\%$ lowest or the $33\%$ highest concentration halos yields $\sigma_{f_{\rm NL}} = 44, 165, 19$, respectively. Our results suggest taking host halo concentration into account in galaxy selection strategies to maximize the signal-to-noise on $f_{\rm NL}$. They also motivate more simulation-based efforts to study the $b_\phi(b_1)$ relation of halos and galaxies., Comment: 21 pages, 4 figures, comments welcome; V2 minor modifications throughout, published in JCAP

The scale-dependent bias effect on the galaxy power spectrum is a very promising probe of the local primordial non-Gaussianity (PNG) parameter $f_{\rm NL}$, but the amplitude of the effect is proportional to $f_{\rm NL}b_{\phi}$, where $b_{\phi}$ is the linear PNG galaxy bias parameter. Our knowledge of $b_{\phi}$ is currently very limited, yet nearly all existing $f_{\rm NL}$ constraints and forecasts assume precise knowledge for it. Here, we use the BOSS DR12 galaxy power spectrum to illustrate how our uncertain knowledge of $b_{\phi}$ currently prevents us from constraining $f_{\rm NL}$ with a given statistical precision $\sigma_{f_{\rm NL}}$. Assuming different fixed choices for the relation between $b_{\phi}$ and the linear density bias $b_1$, we find that $\sigma_{f_{\rm NL}}$ can vary by as much as an order of magnitude. Our strongest bound is $f_{\rm NL} = 16 \pm 16\ (1\sigma)$, while the loosest is $f_{\rm NL} = 230 \pm 226\ (1\sigma)$ for the same BOSS data. The impact of $b_{\phi}$ can be especially pronounced because it can be close to zero. We also show how marginalizing over $b_{\phi}$ with wide priors is not conservative, and leads in fact to biased constraints through parameter space projection effects. Independently of galaxy bias assumptions, the scale-dependent bias effect can only be used to detect $f_{\rm NL} \neq 0$ by constraining the product $f_{\rm NL}b_{\phi}$, but the error bar $\sigma_{f_{\rm NL}}$ remains undetermined and the results cannot be compared with the CMB; we find $f_{\rm NL}b_{\phi} \neq 0$ with $1.6\sigma$ significance. We also comment on why these issues are important for analyses with the galaxy bispectrum. Our results strongly motivate simulation-based research programs aimed at robust theoretical priors for the $b_{\phi}$ parameter, without which we may never be able to competitively constrain $f_{\rm NL}$ using galaxy data., Comment: 15 pages, 6 figures, 3 tables. Comments welcomed! v2 Matches version published in JCAP

The integrated shear 3-point correlation function $\zeta_{\pm}$ is a higher-order statistic of the cosmic shear field that describes the modulation of the 2-point correlation function $\xi_{\pm}$ by long-wavelength features in the field. Here, we introduce a new theoretical model to calculate $\zeta_{\pm}$ that is accurate on small angular scales, and that allows to take baryonic feedback effects into account. Our model builds on the realization that the small-scale $\zeta_{\pm}$ is dominated by the nonlinear matter bispectrum in the squeezed limit, which can be evaluated accurately using the nonlinear matter power spectrum and its first-order response functions to density and tidal field perturbations. We demonstrate the accuracy of our model by showing that it reproduces the small-scale $\zeta_{\pm}$ measured in simulated cosmic shear maps. The impact of baryonic feedback enters effectively only through the corresponding impact on the nonlinear matter power spectrum, thereby permitting to account for these astrophysical effects on $\zeta_{\pm}$ similarly to how they are currently accounted for on $\xi_{\pm}$. Using a simple idealized Fisher matrix forecast for a DES-like survey we find that, compared to $\xi_{\pm}$, a combined $\xi_{\pm}\ \&\ \zeta_{\pm}$ analysis can lead to improvements of order $20-40\%$ on the constraints of cosmological parameters such as $\sigma_8$ or the dark energy equation of state parameter $w_0$. We find similar levels of improvement on the constraints of the baryonic feedback parameters, which strengthens the prospects for cosmic shear data to obtain tight constraints not only on cosmology but also on astrophysical feedback models. These are encouraging results that motivate future works on the integrated shear 3-point correlation function towards applications to real survey data., Comment: 15 pages, 9 figures; v2 matches the published version in MNRAS

We use separate universe simulations with the IllustrisTNG galaxy formation model to predict the local PNG bias parameters $b_\phi$ and $b_{\phi\delta}$ of atomic neutral hydrogen, ${\rm H_I}$. These parameters and their relation to the linear density bias parameter $b_1$ play a key role in observational constraints of the local PNG parameter $f_{\rm NL}$ using the ${\rm H_I}$ power spectrum and bispectrum. Our results show that the popular calculation based on the universality of the halo mass function overpredicts the $b_\phi(b_1)$ and $b_{\phi\delta}(b_1)$ relations measured in the simulations. In particular, our results show that at $z \lesssim 1$ the ${\rm H_I}$ power spectrum is more sensitive to $f_{\rm NL}$ compared to previously thought ($b_\phi$ is more negative), but is less sensitive at other epochs ($b_\phi$ is less positive). We discuss how this can be explained by the competition of physical effects such as that large-scale gravitational potentials with local PNG (i) accelerate the conversion of hydrogen to heavy elements by star formation, (ii) enhance the effects of baryonic feedback that eject the gas to regions more exposed to ionizing radiation, and (iii) promote the formation of denser structures that shield the ${\rm H_I}$ more efficiently. Our numerical results can be used to revise existing forecast studies on $f_{\rm NL}$ using 21cm line-intensity mapping data. Despite this first step towards predictions for the local PNG bias parameters of ${\rm H_I}$, we emphasize that more work is needed to assess their sensitivity on the assumed galaxy formation physics and ${\rm H_I}$ modeling strategy., Comment: 13 pages, 5 figures, 1 table. Comments welcomed! v2 matches version published in JCAP

We use hydrodynamical separate universe simulations with the IllustrisTNG model to predict the local primordial non-Gaussianity (PNG) bias parameters $b_{\phi}$ and $b_{\phi\delta}$, which enter at leading order in the galaxy power spectrum and bispectrum. This is the first time that $b_{\phi\delta}$ is measured from either gravity-only or galaxy formation simulations. For dark matter halos, the popular assumption of universality overpredicts the $b_{\phi\delta}(b_1)$ relation in the range $1 \lesssim b_1 \lesssim 3$ by up to $\Delta b_{\phi\delta} \sim 3$ ($b_1$ is the linear density bias). The adequacy of the universality relation is worse for the simulated galaxies, with the relations $b_{\phi}(b_1)$ and $b_{\phi\delta}(b_1)$ being generically redshift-dependent and very sensitive to how galaxies are selected (we test total, stellar and black hole mass, black hole mass accretion rate and color). The uncertainties on $b_{\phi}$ and $b_{\phi\delta}$ have a direct, often overlooked impact on the constraints of the local PNG parameter $f_{\rm NL}$, which we study and discuss. For a survey with $V = 100{\rm Gpc}^3/h^3$ at $z=1$, uncertainties $\Delta b_{\phi} \lesssim 1$ and $\Delta b_{\phi\delta} \lesssim 5$ around values close to the fiducial can yield relatively unbiased constraints on $f_{\rm NL}$ using power spectrum and bispectrum data. We also show why priors on galaxy bias are useful even in analyses that fit for products $f_{\rm NL} b_{\phi}$ and $f_{\rm NL} b_{\phi\delta}$. The strategies we discuss to deal with galaxy bias uncertainties can be straightforwardly implemented in existing $f_{\rm NL}$ constraint analyses (we provide fits for some of the bias relations). Our results motivate more works with galaxy formation simulations to refine our understanding of $b_{\phi}$ and $b_{\phi\delta}$ towards improved constraints on $f_{\rm NL}$., Comment: 18 pages + 2 appendices; 8 figures. Comments welcomed! v2: added references and the new Fig. 6 with an additional forecast (same conclusions). v3: minor comments added following the referee report

We use the forward modeling approach to galaxy clustering combined with the likelihood from the effective-field theory of large-scale structure to measure assembly bias, i.e. the dependence of halo bias on properties beyond the total mass, in the linear ($b_1$) and second order bias parameters ($b_2$ and $b_{K^2}$) of dark matter halos in $N$-body simulations. This is the first time that assembly bias in the tidal bias parameter $b_{K^2}$ is measured. We focus on three standard halo properties: the concentration $c$, spin $\lambda$, and sphericity $s$, for which we find an assembly bias signal in $b_{K^2}$ that is opposite to that in $b_1$. Specifically, at fixed mass, halos that get more (less) positively biased in $b_1$, get less (more) negatively biased in $b_{K^2}$. We also investigate the impact of assembly bias on the $b_2(b_1)$ and $b_{K^2}(b_1)$ relations, and find that while the $b_2(b_1)$ relation stays roughly unchanged, assembly bias strongly impacts the $b_{K^2}(b_1)$ relation. This impact likely extends also to the corresponding relation for galaxies, which motivates future studies to design better priors on $b_{K^2}(b_1)$ for use in cosmological constraints from galaxy clustering data., Comment: 22 pages, 5 figures ; v2: version published in JCAP

We use field-level forward models of galaxy clustering and the EFT likelihood formalism to study, for the first time for self-consistently simulated galaxies, the relations between the linear $b_1$ and second-order bias parameters $b_2$ and $b_{K^2}$. The forward models utilize all of the information available in the galaxy distribution up to a given order in perturbation theory, which allows us to infer these bias parameters with high signal-to-noise, even from relatively small volumes ($L_{\rm box} = 205{\rm Mpc}/h$). We consider galaxies from the IllustrisTNG simulations, and our main result is that the $b_2(b_1)$ and $b_{K^2}(b_1)$ relations obtained from gravity-only simulations for total mass selected objects are broadly preserved for simulated galaxies selected by stellar mass, star formation rate, color and black hole accretion rate. We also find good agreement between the bias relations of the simulated galaxies and a number of recent estimates for observed galaxy samples. The consistency under different galaxy selection criteria suggests that theoretical priors on these bias relations may be used to improve cosmological constraints based on observed galaxy samples. We do identify some small differences between the bias relations in the hydrodynamical and gravity-only simulations, which we show can be linked to the environmental dependence of the relation between galaxy properties and mass. We also show that the EFT likelihood recovers the value of $\sigma_8$ to percent-level from various galaxy samples (including splits by color and star formation rate) and after marginalizing over 8 bias parameters. This demonstration using simulated galaxies adds to previous works based on halos as tracers, and strengthens further the potential of forward models to infer cosmology from galaxy data., Comment: 19 pages, 9 figures. Version 2 matches published version and contains a comparison to bias estimates from observations. Comments welcomed!

Halo occupation distribution (HOD) models describe the number of galaxies that reside in different haloes, and are widely used in galaxy-halo connection studies using the halo model (HM). Here, we introduce and study HOD response functions $R_\mathcal{O}^g$ that describe the response of the HODs to long-wavelength perturbations $\mathcal{O}$. The linear galaxy bias parameters $b_\mathcal{O}^g$ are a weighted version of $b_\mathcal{O}^h + R_\mathcal{O}^g$, where $b_\mathcal{O}^h$ is the halo bias, but the contribution from $R_\mathcal{O}^g$ is routinely ignored in the literature. We investigate the impact of this by measuring the $R_\mathcal{O}^g$ in separate universe simulations of the IllustrisTNG model for three types of perturbations: total matter perturbations, $\mathcal{O}=\delta_m$; baryon-CDM compensated isocurvature perturbations, $\mathcal{O}=\sigma$; and potential perturbations with local primordial non-Gaussianity, $\mathcal{O}\propto f_{\rm NL}\phi$. Our main takeaway message is that the $R_\mathcal{O}^g$ are not negligible in general and their size should be estimated on a case-by-case basis. For stellar-mass selected galaxies, the responses $R_\phi^g$ and $R_\sigma^g$ are sizeable and cannot be neglected in HM calculations of the bias parameters $b_\phi^g$ and $b_\sigma^g$; this is relevant to constrain inflation using galaxies. On the other hand, we do not detect a strong impact of the HOD response $R_1^g$ on the linear galaxy bias $b_1^g$. These results can be explained by the impact that the perturbations have on stellar-to-total-mass relations. We also look into the impact on the bias of the gas distribution and find similar conclusions. We show that a single extra parameter describing the overall amplitude of $R_\mathcal{O}^g$ recovers the measured $b_\mathcal{O}^g$ well, which indicates that $R_\mathcal{O}^g$ can be easily added to HM/HOD studies as a new ingredient., Comment: 22 pages, 12 figures, 1 table. Comments are welcome! Accepted by JCAP

We study the impact that uncertainties on assumed relations between galaxy bias parameters have on constraints of the local PNG $f_{\rm NL}$ parameter. We focus on the relation between the linear density galaxy bias $b_1$ and local PNG bias $b_\phi$ in an idealized forecast setup with multitracer galaxy power spectrum and bispectrum data. We consider two parametrizations of galaxy bias: 1) one inspired by the universality relation where $b_\phi = 2\delta_c\left(b_1 - p\right)$ and $p$ is a free parameter; and 2) another in which the product of bias parameters and $f_{\rm NL}$, like $f_{\rm NL} b_\phi$, is directly fitted for. The constraints on the $f_{\rm NL}-p$ plane are markedly bimodal, and both the central value and width of marginalized constraints on $f_{\rm NL}$ depend sensitively on the priors on $p$. Assuming fixed $p=1$ in the constraints with a fiducial value of $p=0.55$ can bias the inferred $f_{\rm NL}$ by $0.5\sigma$ to $1\sigma$; priors $\Delta p \approx 0.5$ around this fiducial value are however sufficient in our setup to return unbiased constraints. In power spectrum analyses, parametrization 2, that makes no assumptions on $b_\phi$, can distinguish $f_{\rm NL} \neq 0$ with the same significance as parametrization 1 assuming perfect knowledge of $b_\phi$ (the value of $f_{\rm NL}$ is however left unknown). A drawback of parametrization 2 is that the addition of the bispectrum information is not as beneficial as in parametrization 1. Our results motivate strongly the incorporation of mitigation strategies for bias uncertainties in PNG constraint analyses, as well as further theoretical studies on the relations between bias parameters to better inform those strategies., Comment: 16 pages + 2 appendices; 7 figures; comments welcomed!

We study the impact that large-scale perturbations of (i) the matter density and (ii) the primordial gravitational potential with local primordial non-Gaussianity (PNG) have on galaxy formation using the IllustrisTNG model. We focus on the linear galaxy bias $b_1$ and the coefficient $b_\phi$ of the scale-dependent bias induced by PNG, which describe the response of galaxy number counts to these two types of perturbations, respectively. We perform our study using separate universe simulations, in which the effect of the perturbations is mimicked by changes to the cosmological parameters: modified cosmic matter density for $b_1$ and modified amplitude $\mathcal{A}_s$ of the primordial scalar power spectrum for $b_\phi$. We find that the widely used universality relation $b_\phi = 2\delta_c(b_1 - 1)$ is a poor description of the bias of haloes and galaxies selected by stellar mass $M_*$, which is instead described better by $b_\phi(M_*) = 2\delta_c(b_1(M_*) - p)$ with $p \in [0.4, 0.7]$. This is explained by the different impact that matter overdensities and local PNG have on the median stellar-to-halo-mass relation. A simple model of this impact allows us to describe the stellar mass dependence of $b_1$ and $b_\phi$ fairly well. Our results also show a nontrivial relation between $b_1$ and $b_\phi$ for galaxies selected by color and black hole mass accretion rate. Our results provide refined priors on $b_\phi$ for local PNG constraints and forecasts using galaxy clustering. Given that the widely used universality relation underpredicts $b_\phi(M_*)$, existing analyses may underestimate the true constraining power on local PNG., Comment: 21 pages, 12 figures, 1 table. Comments welcomed! v2: version published in JCAP

We investigate the potential of the galaxy power spectrum to constrain compensated isocurvature perturbations (CIPs), primordial fluctuations in the baryon density that are compensated by fluctuations in CDM density to ensure an unperturbed total matter density. We show that CIPs contribute to the galaxy overdensity at linear order, and if they are close to scale-invariant, their effects are nearly perfectly degenerate with the local PNG parameter $f_{\rm nl}$ if they correlate with the adiabatic perturbations. This degeneracy can however be broken by analyzing multiple galaxy samples with different bias parameters, or by taking CMB priors on $f_{\rm nl}$ into account. Parametrizing the amplitude of the CIP power spectrum as $P_{\sigma\sigma} = A^2P_{\mathcal{R}\mathcal{R}}$ (where $P_{\mathcal{R}\mathcal{R}}$ is the adiabatic power spectrum) we find, for a number of fiducial galaxy samples in a simplified forecast setup, that constraints on $A$, relative to those on $f_{\rm nl}$, of order $\sigma_{A}/\sigma_{f_{\rm nl}} \approx 1-2$ are achievable for CIPs correlated with adiabatic perturbations, and $\sigma_{A}/\sigma_{f_{\rm nl}} \approx 5$ for the uncorrelated case. These values are independent of survey volume, and suggest that current galaxy data are already able to improve significantly on the tightest existing constraints on CIPs from the CMB. Future galaxy surveys that aim to achieve $\sigma_{f_{\rm nl}} \sim 1$ have the potential to place even stronger bounds on CIPs., Comment: 17 pages + references and appendix; 6 figures; 1 table. Comments welcomed! v2 matches version published in JCAP

We investigate the accuracy of weak lensing simulations by comparing the results of five independently developed lensing simulation codes run on the same input $N$-body simulation. Our comparison focuses on the lensing convergence maps produced by the codes, and in particular on the corresponding PDFs, power spectra and peak counts. We find that the convergence power spectra of the lensing codes agree to $\lesssim 2\%$ out to scales $\ell \approx 4000$. For lensing peak counts, the agreement is better than $5\%$ for peaks with signal-to-noise $\lesssim 6$. We also discuss the systematic errors due to the Born approximation, line-of-sight discretization, particle noise and smoothing. The lensing codes tested deal in markedly different ways with these effects, but they nonetheless display a satisfactory level of agreement. Our results thus suggest that systematic errors due to the operation of existing lensing codes should be small. Moreover their impact on the convergence power spectra for a lensing simulation can be predicted given its numerical details, which may then serve as a validation test., Comment: 16 pages, 10 figures; comments welcomed! In v2: minor clarifications added; matches version published in MNRAS

The large-scale clustering of matter is impacted by baryonic physics, particularly AGN feedback. Modelling or mitigating this impact will be essential for making full use of upcoming measurements of cosmic shear and other large-scale structure probes. We study baryonic effects on the matter bispectrum, using measurements from a selection of state-of-the-art hydrodynamical simulations: IllustrisTNG, Illustris, EAGLE, and BAHAMAS. We identify a low-redshift enhancement of the bispectrum, peaking at $k\sim 3h\,{\rm Mpc}^{-1}$, that is present in several simulations, and discuss how it can be associated to the evolving nature of AGN feedback at late times. This enhancement does not appear in the matter power spectrum, and therefore represents a new source of degeneracy breaking between two- and three-point statistics. In addition, we provide physical interpretations for other aspects of these measurements, and make initial comparisons to predictions from perturbation theory, empirical fitting formulas, and the response function formalism. We publicly release our measurements (including estimates of their uncertainty due to sample variance) and bispectrum measurement code as resources for the community., Comment: 29 pages, 24 figures. v2: MNRAS published version, with references and clarifications added (conclusions unchanged). See Fig. 13 for a summary of the main results. Simulation measurements are available at https://github.com/sjforeman/hydro_bispectrum , and bispectrum measurement code is available at https://github.com/sjforeman/bskit

We introduce The Novel Probes Project, an initiative to advance the field of astrophysical tests of the dark sector by creating a forum that connects observers and theorists. This review focuses on tests of gravity and is intended to be of use primarily to observers, but also to theorists with interest in the development of experimental tests. It is twinned with a separate review on tests of dark matter self-interactions (Adhikari et al., in prep.). Our focus is on astrophysical probes of gravity in the weak-field regime, ranging from stars to quasilinear cosmological scales. These are complementary to both strong-field tests and background and linear probes in cosmology. In particular, the nonlinear screening mechanisms that are an integral part of viable modified gravity models lead to characteristic signals specifically on astrophysical scales. The constraining power of these signals is not limited by cosmic variance, but comes with the challenge of building robust theoretical models of the nonlinear dynamics of stars, galaxies, clusters and large scale structure. In this review we lay the groundwork for a thorough exploration of the astrophysical regime with an eye to using the current and next generation of observations for tests of gravity. We begin by setting the scene for how theories beyond General Relativity are expected to behave, focusing primarily on screened fifth forces. We describe the analytic and numerical techniques for exploring the pertinent astrophysical systems, as well as the signatures of modified gravity. With these in hand we present a range of observational tests, and discuss prospects for future measurements and theoretical developments., Comment: 108 pages; matches Reviews of Modern Physics accepted version; project website at https://www.novelprobes.org/

We study the impact that baryon-CDM relative density perturbations $\delta_{bc}$ have on galaxy formation using cosmological simulations with the IllustrisTNG model. These isocurvature (non-adiabatic) perturbations can be induced primordially, if multiple fields are present during inflation, and are generated before baryon-photon decoupling when baryons did not comove with CDM. The presence of long-wavelength $\delta_{bc}$ perturbations in our simulations is mimicked by modifying the ratios of the cosmic densities of baryons $\Omega_b$ and CDM $\Omega_c$, at fixed total matter density $\Omega_m$. We measure the corresponding galaxy bias parameter $b_{\delta}^{bc}$ as the response of galaxy abundances to $\delta_{bc}$. When selecting by total host halo mass, $b_{\delta}^{bc}$ is negative and it decreases with mass and redshift. Stellar-mass selected simulated galaxies show a weaker or even the opposite trend because of the competing effects of $\delta_{bc}$ on the halo mass function and stellar-to-halo-mass relations. We show that simple modeling of the latter two effects describes $b_{\delta}^{bc}$ for stellar-mass-selected objects well. We find $b_{\delta}^{bc} =0.6$ for $M_* = 10^{11}\ M_{\odot}/h$ and $z=0.5$, which is representative of BOSS DR12 galaxies. For $\delta_{bc}$ modes generated by baryon-photon interactions, we estimate the impact on the DR12 power spectrum to be below $1\%$, and shifts on inferred distance and growth rate parameters should not exceed $0.1\%$., Comment: 18 pages (+8 pages with appendices and references); 8 figures; 1 table; comments welcomed. In v2: (i) extended discussion with new figure on the generation of baryon-CDM perturbations; (ii) new figure quantifying the contributions to the bias parameter from changes in the mass function and changes in the stellar-to-halo-mass relations. v3 matches version published in JCAP

We measure power spectrum response functions in the presence of baryonic physical processes using separate universe simulations with the IllustrisTNG galaxy formation model. The response functions describe how the small-scale power spectrum reacts to long-wavelength perturbations and they can be efficiently measured with the separate universe technique by absorbing the effects of the long modes into a modified cosmology. Specifically, we focus on the total first-order matter power spectrum response to an isotropic density fluctuation $R_1(k,z)$, which is fully determined by the logarithmic derivative of the nonlinear matter power spectrum ${\rm dln}P_m(k,z)/{\rm dln}k$ and the growth-only response function $G_1(k,z)$. We find that $G_1(k,z)$ is not affected by the baryonic physical processes in the simulations at redshifts $z < 3$ and on all scales probed ($k \lesssim 15h/{\rm Mpc}$, i.e. length scales $\gtrsim 0.4 {\rm Mpc}/h$). In practice, this implies that the power spectrum fully specifies the baryonic dependence of its response function. Assuming an idealized lensing survey setup, we evaluate numerically the baryonic impact on the squeezed-lensing bispectrum and the lensing super-sample power spectrum covariance, which are given in terms of responses. Our results show that these higher-order lensing statistics can display varying levels of sensitivity to baryonic effects compared to the power spectrum, with the squeezed-bispectrum being the least sensitive. We also show that ignoring baryonic effects on lensing covariances slightly overestimates the error budget (and is therefore conservative from the point of view of parameter error bars) and likely has negligible impact on parameter biases in inference analyses., Comment: 15 pages, 6 figures, 1 table; comments welcomed! v2 matches version published in MNRAS

We present a calculation of the angle-averaged squeezed matter bispectrum covariance ${\rm Cov}\left(B_{m}(k_1, k_1', s_1), B_{m}(k_2, k_2', s_2)\right)$, $s_i \ll k_i,k_i'$ ($i=1,2$), that uses matter power spectrum responses to describe the coupling of large- to short-scale modes in the nonlinear regime. The covariance is given by a certain configuration of the 6-point function, which we show is dominated by response-type mode-coupling terms in the squeezed bispectrum limit. The terms that are not captured by responses are small, effectively rendering our calculation complete and predictive for linear $s_1,s_2$ values and any nonlinear values of $k_1,k_1',k_2,k_2'$. Our numerical results show that the squeezed bispectrum super-sample covariance is only a negligible contribution. We also compute the power spectrum-bispectrum cross-covariance using responses. Our derivation for the squeezed matter bispectrum is the starting point to calculate analytical covariances for more realistic galaxy clustering and weak-lensing applications. It can also be used in cross-checks of numerical ensemble estimates of the general bispectrum covariance, given that it is effectively noise-free and complete in the squeezed limit., Comment: 29 pages + appendices; 3 figures and 1 table. The busy reader can focus on the summary in the Introduction, Table 1 and Section 6 with results. Comments welcomed. v2: minor corrections to the text; version published on JCAP

We present a general method to compute the nonlinear matter power spectrum for dark energy and modified gravity scenarios with percent-level accuracy. By adopting the halo model and nonlinear perturbation theory, we predict the reaction of a $\Lambda$CDM matter power spectrum to the physics of an extended cosmological parameter space. By comparing our predictions to $N$-body simulations we demonstrate that with no-free parameters we can recover the nonlinear matter power spectrum for a wide range of different $w_0$-$w_a$ dark energy models to better than 1% accuracy out to $k \approx 1 \, h \, {\rm Mpc}^{-1}$. We obtain a similar performance for both DGP and $f(R)$ gravity, with the nonlinear matter power spectrum predicted to better than 3% accuracy over the same range of scales. When including direct measurements of the halo mass function from the simulations, this accuracy improves to 1%. With a single suite of standard $\Lambda$CDM $N$-body simulations, our methodology provides a direct route to constrain a wide range of non-standard extensions to the concordance cosmology in the high signal-to-noise nonlinear regime., Comment: 18+6 pages, 15 figures. Matches MNRAS accepted version. Few typos corrected

Accurate inference of cosmology from weak lensing shear requires an accurate shear power spectrum covariance matrix. Here, we investigate this accuracy requirement and quantify the relative importance of the Gaussian (G), super-sample covariance (SSC) and connected non-Gaussian (cNG) contributions to the covariance. Specifically, we forecast cosmological parameter constraints for future wide-field surveys and study how different covariance matrix components affect parameter bounds. Our main result is that the cNG term represents only a small and potentially negligible contribution to statistical parameter errors: the errors obtained using the G+SSC subset are within $\lesssim 5\%$ of those obtained with the full G+SSC+cNG matrix for a Euclid-like survey. This result also holds for the shear two-point correlation function, variations in survey specifications and for different analytical prescriptions of the cNG term. The cNG term is that which is often tackled using numerically expensive ensembles of survey realizations. Our results suggest however that the accuracy of analytical or approximate numerical methods to compute the cNG term is likely to be sufficient for cosmic shear inference from the next generation of surveys., Comment: 15 pages, 5 figures, 1 table. Comments welcomed. Changes in v2: (i) improved visualization of Fig.1; (ii) added one more test of varying survey specifications

We derive the complete super-sample covariance (SSC) of the matter and weak lensing convergence power spectra using the power spectrum response formalism to accurately describe the coupling of super- to sub-survey modes. The SSC term is completely characterized by the survey window function, the nonlinear matter power spectrum and the full first-order nonlinear power spectrum response function, which describes the response to super-survey density and tidal field perturbations. Generalized separate universe simulations can efficiently measure these responses in the nonlinear regime of structure formation, which is necessary for lensing applications. We derive the lensing SSC formulae for two cases: one under the Limber and flat-sky approximations, and a more general one that goes beyond the Limber approximation in the super-survey mode and is valid for curved sky applications. Quantitatively, we find that for sky fractions $f_{\rm sky} \approx 0.3$ and a single source redshift at $z_S=1$, the use of the flat-sky and Limber approximation underestimates the total SSC contribution by $\approx 10\%$. The contribution from super-survey tidal fields to the lensing SSC, which has not been included in cosmological analyses so far, is shown to represent about $5\%$ of the total lensing covariance on multipoles $\ell_1,\ell_2 \gtrsim 300$. The SSC is the dominant off-diagonal contribution to the total lensing covariance, making it appropriate to include these tidal terms and beyond flat-sky/Limber corrections in cosmic shear analyses., Comment: 23 pages + appendices; 5 figures. Major improvements in v3: (i) added derivation of SSC beyond flat-sky/Limber approximations; (ii) results now use simulation measurements of the tidal response coefficient. Conclusions remain otherwise the same

Cosmological models with Galileon gravity are an alternative to the standard $\Lambda {\rm CDM}$ paradigm with testable predictions at the level of its self-accelerating solutions for the expansion history, as well as large-scale structure formation. Here, we place constraints on the full parameter space of these models using data from the cosmic microwave background (CMB) (including lensing), baryonic acoustic oscillations (BAO) and the Integrated Sachs-Wolfe (ISW) effect. We pay special attention to the ISW effect for which we use the cross-spectra, $C_\ell^{\rm T g}$, of CMB temperature maps and foreground galaxies from the WISE survey. The sign of $C_\ell^{\rm T g}$ is set by the time evolution of the lensing potential in the redshift range of the galaxy sample: it is positive if the potential decays (like in $\Lambda {\rm CDM}$), negative if it deepens. We constrain three subsets of Galileon gravity separately known as the Cubic, Quartic and Quintic Galileons. The cubic Galileon model predicts a negative $C_\ell^{\rm T g}$ and exhibits a $7.8\sigma$ tension with the data, which effectively rules it out. For the quartic and quintic models the ISW data also rule out a significant region of the parameter space but permit regions where the goodness-of-fit is comparable to $\Lambda {\rm CDM}$. The data prefers a non zero sum of the neutrino masses ($\sum m_\nu\approx 0.5$eV) with $ \sim \! 5\sigma$ significance in these models. The best-fitting models have values of $H_0$ consistent with local determinations, thereby avoiding the tension that exists in $\Lambda {\rm CDM}$. We also identify and discuss a $\sim \! 2\sigma$ tension that Galileon gravity exhibits with recent BAO measurements. Our analysis shows overall that Galileon cosmologies cannot be ruled out by current data but future lensing, BAO and ISW data hold strong potential to do so., Comment: 26 pages, 7 Figures; published version

We present a calculation of the matter power spectrum covariance matrix ${\rm Cov}(\bf{k}_1,\bf{k}_2)$ that uses power spectrum responses to accurately describe the coupling between large- and small-scale modes beyond the perturbative regime. These response functions can be measured with (small-volume) N-body simulations, which is why the response contributions to the covariance remain valid and predictive at all orders in perturbation theory. A novel and key step presented here is the use of responses to compute loop contributions with soft loop momenta, which extends the application of the response approach beyond that of previously considered squeezed $n$-point functions. The calculation presented here does not involve any fitting parameters. When including response-type terms up to 1-loop order in perturbation theory, we find that our calculation captures the bulk of the total covariance as estimated from simulations up to values of $k_1, k_2 \sim 1\ h/{\rm Mpc}$. Moreover, the prediction is guaranteed to be accurate whenever the softer mode is sufficiently linear, ${\rm min}\{k_1,k_2\} \lesssim 0.08\ h/{\rm Mpc}$. We identify and discuss straightforward improvements in the context of the response approach, which are expected to further increase the accuracy of the calculation presented here., Comment: 26 pages + appendices and references. 10 figures. Minor change in caption of Fig.3 in v2

We introduce a rigorous definition of general power-spectrum responses as resummed vertices with two hard and $n$ soft momenta in cosmological perturbation theory. These responses measure the impact of long-wavelength perturbations on the local small-scale power spectrum. The kinematic structure of the responses (i.e., their angular dependence) can be decomposed unambiguously through a "bias" expansion of the local power spectrum, with a fixed number of physical response coefficients, which are only a function of the hard wavenumber $k$. Further, the responses up to $n$-th order completely describe the $(n+2)$-point function in the squeezed limit, i.e. with two hard and $n$ soft modes, which one can use to derive the response coefficients. This generalizes previous results, which relate the angle-averaged squeezed limit to isotropic response coefficients. We derive the complete expression of first- and second-order responses at leading order in perturbation theory, and present extrapolations to nonlinear scales based on simulation measurements of the isotropic response coefficients. As an application, we use these results to predict the non-Gaussian part of the angle-averaged matter power spectrum covariance ${\rm Cov}^{\rm NG}_{\ell = 0}(k_1,k_2)$, in the limit where one of the modes, say $k_2$, is much smaller than the other. Without any free parameters, our model results are in very good agreement with simulations for $k_2 \lesssim 0.06\ h/{\rm Mpc}$, and for any $k_1 \gtrsim 2 k_2$. The well-defined kinematic structure of the power spectrum response also permits a quick evaluation of the angular dependence of the covariance matrix. While we focus on the matter density field, the formalism presented here can be generalized to generic tracers such as galaxies., Comment: 28 pages + appendices and references. 5 figures. Comments welcomed

We describe and demonstrate the potential of a new and very efficient method for simulating certain classes of modified gravity theories, such as the widely studied $f(R)$ gravity models. High resolution simulations for such models are currently very slow due to the highly nonlinear partial differential equation that needs to be solved exactly to predict the modified gravitational force. This nonlinearity is partly inherent, but is also exacerbated by the specific numerical algorithm used, which employs a variable redefinition to prevent numerical instabilities. The standard Newton-Gauss-Seidel iterative method used to tackle this problem has a poor convergence rate. Our new method not only avoids this, but also allows the discretised equation to be written in a form that is analytically solvable. We show that this new method greatly improves the performance and efficiency of $f(R)$ simulations. For example, a test simulation with $512^3$ particles in a box of size $512 \, \mathrm{Mpc}/h$ is now 5 times faster than before, while a Millennium-resolution simulation for $f(R)$ gravity is estimated to be more than 20 times faster than with the old method. Our new implementation will be particularly useful for running very high resolution, large-sized simulations which, to date, are only possible for the standard model, and also makes it feasible to run large numbers of lower resolution simulations for covariance analyses. We hope that the method will bring us to a new era for precision cosmological tests of gravity., Comment: 14 pages, 3 figures, 1 table; version accepted for publication in JCAP

We present high signal-to-noise galaxy-galaxy lensing measurements of the BOSS CMASS sample using 250 square degrees of weak lensing data from CFHTLenS and CS82. We compare this signal with predictions from mock catalogs trained to match observables including the stellar mass function and the projected and two dimensional clustering of CMASS. We show that the clustering of CMASS, together with standard models of the galaxy-halo connection, robustly predicts a lensing signal that is 20-40% larger than observed. Detailed tests show that our results are robust to a variety of systematic effects. Lowering the value of $S_{\rm 8}=\sigma_{\rm 8} \sqrt{\Omega_{\rm m}/0.3}$ compared to Planck2015 reconciles the lensing with clustering. However, given the scale of our measurement ($r<10$ $h^{-1}$ Mpc), other effects may also be at play and need to be taken into consideration. We explore the impact of baryon physics, assembly bias, massive neutrinos, and modifications to general relativity on $\Delta\Sigma$ and show that several of these effects may be non-negligible given the precision of our measurement. Disentangling cosmological effects from the details of the galaxy-halo connection, the effects of baryons, and massive neutrinos, is the next challenge facing joint lensing and clustering analyses. This is especially true in the context of large galaxy samples from Baryon Acoustic Oscillation surveys with precise measurements but complex selection functions., Comment: 26 pages. Submitted to MNRAS. Comments welcome

We study the imprints that theories of gravity beyond GR can leave on the lensing signal around line of sight directions that are predominantly halo-underdense (called troughs) and halo-overdense. To carry out our investigations, we consider the normal branch of DGP gravity, as well as a phenomenological variant thereof that directly modifies the lensing potential. The predictions of these models are obtained with N-body simulation and ray-tracing methods using the ECOSMOG and Ray-Ramses codes. We analyse the stacked lensing convergence profiles around the underdense and overdense lines of sight, which exhibit, respectively, a suppression and a boost w.r.t. the mean in the field of view. The modifications to gravity in these models strengthen the signal w.r.t. $\Lambda{\rm CDM}$ in a scale-independent way. We find that the size of this effect is the same for both underdense and overdense lines of sight, which implies that the density field along the overdense directions on the sky is not sufficiently evolved to trigger the suppression effects of the screening mechanism. These results are robust to variations in the minimum halo mass and redshift ranges used to identify the lines of sight, as well as to different line of sight aperture sizes and criteria for their underdensity and overdensity thresholds., Comment: 20 pages, 17 figures. v2 matches version published in JCAP

We perform a validation of estimates of the growth rate of structure, described by the parameter combination $f\sigma_8$, in modified gravity cosmologies. We consider an analysis pipeline based on the redshift-space distortion modelling of the clustering wedges statistic of the galaxy correlation function and apply it to mock catalogues of $\Lambda{\rm CDM}$ and the normal branch of DGP cosmologies. We employ a halo occupation distribution approach to construct our mocks, which we ensure resemble the CMASS sample from BOSS in terms of the total galaxy number density and large scale amplitude of the power spectrum monopole. We show that the clustering wedges model successfully recovers the true growth rate difference between DGP and $\Lambda{\rm CDM}$, even for cases with over 40\% enhancement in $f\sigma_8$ compared to $\Lambda{\rm CDM}$. The unbiased performance of the clustering wedges model allows us to use the growth rate values estimated from the BOSS DR12 data to constrain the cross-over scale $r_c$ of DGP gravity to $\left[r_cH_0\right]^{-1} < 0.97$ ($2\sigma$) or $r_c > 3090\ {\rm Mpc}/h$, cutting into the interesting region of parameter space with $r_c \sim H_0^{-1}$ using constraints from the growth of structure alone., Comment: 18 pages, 8 figures. Matches published version in PRD that validates the DR12 RSD clustering wedges model. Versions 1 and 2 contain the validation analysis of the DR11 model